Constant flux relation for driven dissipative systems 
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Conservation laws constrain the stationary state statistics of driven dissipative systems because 
the average flux of a conserved quantity between driving and dissipation scales should be constant. 
This requirement leads to a universal scaling law for flux-measuring correlation functions, which 
generalizes the 4/5-th law of Navier-Stokes turbulence. We demonstrate the utility of this simple idea 
by deriving new exact scaling relations for models of aggregating particle systems in the fluctuation- 
dominated regime and for energy and wave action cascades in models of strong wave turbulence. 
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A major challenge of theoretical physics is to develop 
a general formalism for understanding properties of sys- 
tems far from equilibrium. In the absence of such a for- 
malism, exact results play an important role. They act as 
precursors for more general theories and as checkpoints 
for phenomenological theories and numerical methods. 
In this paper we derive an exact relation satisfied by a 
particular correlation function for an important class of 
turbulent non-equilibrium systems. 

By turbulent systems, we understand not just fluid 
turbulence, but a more general class of driven dissipa- 
tive systems characterised by the presence of an inertial 
range. This is a wide range of scales separating the driv- 
ing and dissipation. In the inertial range, dynamics is 
expected to be universal, i. e. independent of details 
of driving and dissipation. In the steady state, the in- 
ertial range dynamics is characterised by constant fluxes 
of conserved quantities from driving to dissipation scale. 
A well known example is Navier-Stokes (NS) turbulence 
Energy injected into the system at a large length 
scale L is dissipated at a length scale I ~ j/ 3 / 4 , where v is 
the kinematic viscosity. In the inertial range ! < r < I, 
there is a constant energy flux directed to small scales. 

If the inertial range dynamics is scale invariant, the 
constant flux of a conserved quantity always fixes the 
scaling of a certain correlation function. An example is 
Kolmogorov's 4/5 law of 3 dimensional NS turbulence 
which states that, in the inertial range, 

<[«!(»*;*) -t>i($*)] 8 > = -\ tr ' (i) 

where vi (r , t) is the longitudinal component of the ve- 
locity at r at time t and e is the average energy flux. 
Pedagogical derivations of this result are in [l], Q . Other 
examples are the Monin-Yaglom relation for passive ad- 
vection (see [H and references therein) and Burgers tur- 
bulence where integrals of motion can be used to prove 
extreme anomalous scaling of all structure functions Q. 

A general law relating scaling properties of correlation 
functions to inertial range conservation laws is still miss- 



ing. As a result, no counterpart of the 4/5-law is known 
for such important turbulent systems such as 3— and 4- 
wave turbulence and cluster-cluster aggregation. In this 
paper we derive a general constant flux relation (CFR) 
which should be obeyed by all flux-determined correla- 
tions of turbulent systems. We start with a heuristic 
argument, followed by an analytic application of this ar- 
gument to the following specific models: 

• Mass model (MM): Consider particles of positive 
mass on a d-dimensional lattice. They evolve in 
time by diffusion, at rate D, aggregation of particles 
at the same site, at rate A(mi,m2), and injection 
of particles of mass too at rate J /mo. For more 
detail, see p. 

• Charge model (CM): The model is similar to MM 
except for the following modifications. The masses 
(now charges) can be positive or negative. Particles 
of charge ±mo & r e input at the same rate J/toq. 

• Wave turbulence (WT): An ensemble of interacting 
dispersive waves maintained far from equilibrium 
by interaction with external sources and sinks. For 
details and a review see @, Q . 

Together, these models relate to a broad range of disci- 
plines, far beyond the common applications of fluid tur- 
bulence. Applications of the MM or CM include sub- 
monolayer epitaxial thin film growth [9(, river networks 
dSEH) force fluctuations in granular bead packs [13] and 
non-equilibrium phase transitions (l3l . [l4[ ■ The constant 
kernel MM maps [IH onto the directed abelian sandpile, 
one of the first models generating power laws from sim- 
ple dynamical rules [l6j |. WT is realized in such diverse 
situations as capillary and gravity waves in liquids, elec- 
tromagnetic waves in non-linear optics and Alfven waves 
in plasmas (see 0, [ItJ for more applications) . 

We present a heuristic argument for a general turbu- 
lent system. It is necessarily abstract, hiding much of the 
physics of specific systems. We believe that the general 
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approach only acquires clarity after application to con- 
crete examples. The idea is very simple. First identify 
the conserved quantity expected to cascade, we shall call 
it /, and the space in which it flows. Second, use the 
equation of motion to write down a Boltzmann-like con- 
tinuity equation for the average density oil in this space. 
This equation identifies a flux-carrying correlation func- 
tion, II, and a nonlinear coupling, T, controlling the re- 
distribution of /. Third, use dimensional analysis to find 
the scaling of II corresponding to constant flux. Suppose 
the flux of / occurs in a <i-dimensional space, spanned 
by p. The degrees of freedom are at least partially in- 
dexed by p. We denote them by up. For example, in 
three dimensional NS turbulence, energy transfer occurs 
between all modes in Fourier space, so d — 3. Also p = k, 
the wave number, and up is v(k), the velocity field. For 
surface waves, motion occurs in three dimensions but en- 
ergy is transferred between wave trains indexed by a two 
dimensional wave-vector, so d = 2. For particle aggre- 
gation problems, mass is transferred in one dimensional 
mass space so d = 1 even though the particles may be 
diffusing on a lattice of arbitrary spatial dimension. 

The system must have an inertial range in p-space, 
within which the dynamics is entirely dominated by the 
/-conserving interactions. Denote the flux by J(p). As- 
sume for simplicity, that the inertial range dynamics is 
isotropic, so we can average over angles in p-space. Let 
pi be the density, and J{p) the flux, of I in p-space, 
where p = |p| . In the inertial range the flux of / through 
a sphere of radius p does not depend on p. Therefore, 
J(p) = C p~ d p. On the other hand, from the equation of 
motion, one can always derive a continuity equation for 
/ in p-space taking the general form: 

/n— 1 
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T is a non-linear coupling and II is a correlation function 
of the up, whose precise forms must be calculated for 
each specific example. In Eq. ^ it is assumed that the 
non-linearity is of degree n — 1 in the up. We take T to be 
homogeneous of degree (3: T(Ap) = A^T(p). Equating 
scaling exponents, we find 

n(p)~ P - nd -P. (3) 

How does this relate to Kolmogorov's 4/5 law? Eq. ([I]) 
is usually derived in real space although the flux is in 
/c-space, so the correspondence is not clear. A fc-space 
derivation of the 4/5 law, along exactly the lines outlined 
here is indeed possible. It can be found in [l8[. Even 
earlier, a similar approach was taken by Kraichnan in his 
seminal work on 2D turbulence [191 ]. To summarise, for 
NS, n = 3 (nonlinearity is quadratic), (3 = 1 (nonlinear- 
ity contains one derivative) and II ~ {vi{ki)vi{k2)vi{k^)) . 
By Eq. ©, LT ~ k~ 1 - 3d . Thus ((v t {r) ~ vi{0)) 3 ) ~ 
r i+3d-3d _ r ^ cons i s tent with Eq. (fTJ) . 



The remainder of this letter is dedicated to the illus- 
tration and analytical verification Eq. ([3]) for the specific 
models listed above. We illustrate the practical meaning 
of the steps outlined above by explicitly identifying II and 
T and, by exact computations, show that Eq. ^ is more 
than dimensional analysis. The examples extend the ap- 
plicability of the 4/5-law beyond the usual cases. MM 
is a pedagogical example of flux conservation outside of 
hydrodynamics. CM illustrates that the idea is also rel- 
evant to statistical conservation laws. WT demonstrates 
the approach for non-quadratic invariants. 

First, consider the MM. The conserved quantity is 
mass. The mass flux is in mass space. Let the reaction 
rate \ mi ,m 2 be a homogeneous function of degree j3. Let 
N(m, x, t) be the local mass distribution. The continu- 
ity equation for the average mass density, m{N(m, x, t)) 
takes the following form in the limit mo — > 0, t — > oo: 

/>oc 

to / dm 1 dm 2 [io ; i,2 - ha.o ~ h-,o,2] = 0, (4) 
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where Iq-,1,2 = A mi , m2 C(mi,m2)£(m — mi - m 2 ) and 
C(mi,m 2 ,t) = (N(mi,S,t)N(m 2 ,x,t) - -^zS( m i ~ 
m2)N{rni,x,t)) is the average number of pairs of par- 
ticles of masses mi and m 2 per lattice site. Ax is the 
lattice spacing. Eq. (j4|) is not supposed to be obvious. A 
detailed derivation for the case j3 — can be found in 0] . 
Let us look for homogeneous solutions to Eq. (j4|) of degree 
h, i.e., C(Ami, Am 2 ) = A ft C(mi, m 2 ). The value of h is 
determined by the balance between the 'plus' and 'mi- 
nus' terms in the left hand side of Eq. ([4]). This balance 
can be made explicit by applying the following Zakharov 
Transformations 0, l20l ] to the 'minus' terms: 



Second integral: 
Third integral: 




after which we obtain 

/>oo />oo 

= / / dm 1 dm 2 I m - mi .m 2 [m y - m\ - mf] (7) 
Jo Jo 

where y = —h — (3 — 3. Eq. ([7]) is satisfied only if y = 1: 
the function /(mi, m 2 ) = (mi + m 2 ) y — m\ ~ m\ is sign 
definite for any y ^ 1 whereas J m , mi , m! is non-negative. 
Therefore the integral in the right hand side of Eq. ([7]) 
is zero only for y — 1. Consequently, h = —3 — (3, is the 
unique scaling solution, under an extra assumption of 
convergence of collision integrals, which will be discussed 
elsewhere. For (3 = 0, the scaling h = —3 — (3 is already 
established. It is trivially true for d > 2 where mean 
field holds, in d = 2 due to cancellation of logarithmic 
corrections @ and in d < 2 by an exact solution [2x1 ] . 

The scaling of the flux correlation function II = 
mCS(m) ~ m -3- ' 3 agrees with general CFR, Eq. 
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for d = 1, n = 3. Note that the result, remarkably, does 
not depend on the transport properties of the particles 
such as mass dependence of the diffusion coefficient. 

Next, consider CM. There are two conserved quanti- 
ties: average charge and average squared charge. Charge 
is conserved in elementary collisions, but squared charge 
is conserved only on average. To see this, consider 
the coagulations between the following pairs of charges: 
(mi,m 2 ), (-TOi,m 2 ), (mi,-m 2 ) and (—mi,— mz). 
They occur with equal probability and the total squared 
charge in the initial and final states are equal. There- 
fore, squared charge is conserved on average. Total flux 
of charge is zero by the symmetry of input. Therefore, 
for the CM, CFR is associated with a flux of squared 
charge. The corresponding continuity equation is 
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where Jq.j 2 = K mi ^ m2 C(mi,m 2 )S[m ± mi ± m 2 ]. Here 
C(mi,m 2 ) has the same meaning as in the MM. Due 
to symmetry C(— mi,m 2 ) = C(mi,m 2 ). We assume 
K mi ,m 2 to be homogeneous of degree (3. Also, S(x ± 
V ±*) = £™,»=o S(x - (-l) m y - (-l) n z). The following 
ZT is applied to Eq. ©: 



. I m mm 2 , 

First integral: (mi,m 2 ) — ► , , (9) 

\mi mi 

/ ffi^ mvri\ 

Second integral: (mi,m 2 ) — > , I , (10) 

\ m 2 m 2 ' 

and integration over m\ is performed to obtain 



k being the wave vector. The equations of motion are 

5U 



w(fc)ag 



(13) 



where U = J dku(k). The equation of motion conserves 
H. There may be additional conservation laws depend- 
ing on the form of u{k). The system is maintained far 
from equilibrium by adding appropriate forcing and dis- 
sipation terms to Eq. (fT3|) . Most commonly, u(k) has 
degree 3 or 4 in Ajj. 

We first study the energy cascade which conserves en- 
ergy in the inertial range. As a result, the average flux of 
energy is constant, which leads to the continuity equation 
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where U is the interaction part of the Hamiltonian. We 
now consider two cases: 3-wave turbulence (n=3) where 

f 2 _ _ _ _ 
u = lW_dk i 8(k-ki-k 2 )T k . klM [aj:a^ i a^ + c.c], (15) 

•* i=i 

and 4-wave turbulence (n=4) where 

u = j Y\dki5{k + ki - k 2 ~ k z )T kM . MM 
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(16) 



In both cases, we assume T to be homogeneous of degree 
(3. For 3-wave turbulence Eq. (fl"4|) . reduces to 
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where y = —h — (3 — 2 and h is the homogeneity exponent 
of C(mi,m 2 ) to be determined. If y = 2, then Eq. (fTTj) 
is satisfied as the integrand vanishes identically due to 
square charge conservation. Generalizing the argument 
made in case of MM, it is possible to check that y = 2 
is the unique scaling solution, given locality. Thus, h = 
—4 — (3. The scaling of the flux measuring correlation 
function II = m 2 C(m)S(m) ~ m -3- ' 3 is consistent with 
Eq. © for d = 1 and n = 3. 

Finally we consider wave turbulence. It is a Hamilto- 
nian system with canonical variables {ag, <%}g gR d: 



H 



dk[ u){k)arat + u(k) 



(12) 



w is a homogeneous function of k and u(k) is the nonlin- 
ear part of the energy density. When u — 0, the Hamilto- 
nian describes free waves with dispersion law u> = co(k), 



X\{dhk 
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[7fc;fci,fc 2 Il0;l,2 — Tk 1 -k,k 2 ^U0,2] — 0' 



(17) 

where no;i, 2 = / Yii=o ^«(^ e ( a fc°fe 1 a fc 2 )) ^ s the fl ux cor- 
relation function. Note an essential difference between 
the flux correlation function for total energy and that 
for quadratic energy. Total energy flux depends on cor- 
relations between fields and time derivatives of fields 
whereas quadratic energy flux would depend on corre- 
lations between fields only. Transforming (ki,k 2 ) — » 
(fc 2 /fci, kk 2 jk\) in the second integral, we obtain 



= J Hidkikf-^Tk^Ilo^ [k y - kf) , 



(18) 



where y = —3d—h — (3 and h is the degree of homogeneity 
of n. Therefore, y = 0, or h = -3d - (3 solves Eq.JTHJ). 
The result n ~ k^ 3 ^ 13 agrees with Eq. ([3]), as the order 
of nonlinearity is n — 3, dimensionality of flux space 
is d and homogeneity index of interaction kernel T is 
(3. It is also consistent with Kolmogorov-Zakharov (KZ) 
theory describing turbulence of weakly non-linear waves 
Q: in this limit the wave action density corresponding 
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to an energy cascade is n (k) ~ k~^ +d \ Within KZ 
theory, II ~ Lulm(aa 2 } ~ uTS(uj)S d {k)n(k) 2 ~ fc- 3rf -' 3 . 
We stress however that, unlike KZ theory, this result is 
equally valid for strong wave turbulence. 

A similar analysis for the 4 wave case yields h = 
—Ad — (3 for the homogeneity degree of the correlation 

function n ,i;2,3 = / IlLo d ^i( Re ( iL k a k 1 H 2 " 1 k 3 ))- Tnis 
is also consistent with KZ theory for weak non-linearity 
and with Eq. ([3]) for the order of non-linearity n = 4. 

In addition to energy, 4-wave turbulence also con- 
serves wave action N = J dkn(k), where (a^a^) = 

n(k\)8{ki — k-i). There is a CFR associated with the 
cascade of wave action. The flux correlation function 

is n ,i ; 2,3 = JnLo^^ TO ( S fc S fe 1 a fc 2 a fc 3 ))- Assume that 
II is a scaling function of degree h. Applying ZT to 
the continuity equation we find This is consistent with 
Eq. (|3|), for order of non-linearity n = 4, dimensional- 
ity of flux space d and homogeneity degree of interaction 
kernel (3. In the limit of small non-linearity, the scal- 
ing h = —Ad — (3 can be also derived from KZ theory: 
for weakly non-linear waves, wave action density in the 
wave action cascade is n(k) ~ ^-(2/3-a+3d)/3^ ^gj-g a 
is dispersion exponent. On the other hand, within KZ 
theory, II - TS d {k)5{u)n(k) 3 ~ fc- 4 ^. 

This ends our illustration of specific examples where 
the CFR idea, expressed heuristically in Eq. ((3]), is used 
to derive new scaling laws respected even by strongly in- 
teracting systems. Probably there are other examples 
waiting to be analysed using the ideas presented in this 
letter. In closing we should emphasise that, although the 
results for specific systems analysed here are exact, the 
CFR is generally not a theorem. Application of the ZT 
implicitly assumes convergence of the collision integral 
on the scaling solution for the flux carrying correlation 
function (in the wave turbulence literature, this is re- 
ferred to as the condition of locality). This assumption 
is usually not readily verifiable, even a posteriori, and will 
require considerable additional effort both theoretically 
and numerically. Nevertheless, there are specific cases 
where these issues are simplified which we would like to 
mention. For the MM, the constant kernel case (/? = 0), 
can be reduced to a differential equation which can be 
explicitly solved. For wave turbulence with certain types 
of ultra-local interaction coefficients, CFR can be derived 
as a theorem by solving some simple recursion relations. 
It is not our intention to present more detailed analysis 
of these particulars here but rather to give an assurance 
that the assumption of locality has a basis in reality. 

To summarize, we demonstrated, using a number of ex- 
amples, that the scaling of a particular correlation func- 
tion can be determined when there is a constant flux of a 
physical quantity in scale-invariant turbulent dynamics. 



This provides an extension of Kolmogorov's 1941 work on 
NS turbulence and Yaglom's work on passive advection 
to a diverse class of systems. 
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